Parametric four-wave mixing toolbox for superconducting resonators 
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We study a superconducting circuit that can act as a toolbox to generate various Bogoliubov-linear 
and nonlinear quantum operations on the microwave photon modes of superconducting resonators 
within one single circuit. The quantum operations are generated by exploring dispersive four- wave 
mixing (FWM) processes involving the resonator modes. Different FWM geometries can be realized 
by adjusting the circuit parameters and by applying appropriate microwave drivings. We illustrate 
this scheme using a circuit made of two superconducting qubits that couple with each other. Each 
qubit couples with one superconducting resonator. We also discuss main sources of quantum errors 
in this system and study the fidelity of the quantum operations by numerical simulation. Our 
scheme provides a practical approach to realize quantum information protocols on superconducting 
resonators. 

PACS numbers: 85.25.Cp, 03.67.Lx, 84.40.Dc 



I. INTRODUCTION 

Superconducting quantum circuits have been inten- 
sively studied as building blocks for quantum informa- 
tion processing^— Among these devices, superconduct- 
ing microwave resonators have demonstrated relatively 
high quality factors and strong coupling with various 
superconducting qubits^ Various quantum optical phe- 
nomena and novel quantum many-body effects involving 
microwave photons, in particular, circuit quantum elec- 
trodynamics (CQED), have been observed in the super- 
conducting resonators. 15 

The generation of non-classical states such as Fock 
states, entangled states, and NOON states in the super- 
conducting resonators has been studied in recent exper- 
iments and theoretical proposals.— ~— In most of these 
schemes, quantum state manipulation is achieved via 
the coupling between superconducting resonators and 
superconducting qubits. Several schemes have been 
proposed to generate unitary transformations on the 
microwave photons by engineering effective interaction 
Hamiltonians between the resonator modes, including 
Kerr and cross-Kerr interactions, beam-splitter opera- 
tion, and squeezing operation.— ~—. In these schemes, the 
proposed circuits can only realize specific quantum oper- 
ation on the resonator modes. In a recent work, Langford 
and coauthors studied an interesting idea of generating 
quantum operations starting from a nonlinear-Kerr in- 
teraction between different optical (or microwave) modes. 
By applying selected classical pumps, full set of quantum 
operations can be generated between selected modest 5 - 

All of the quantum operations that are crucial for im- 
plementing quantum information protocols on resonator 
modes can be constructed from a basic set of quantum op- 
erations. For both discrete-state and continuous variable 
quantum information protocols on the resonator modes, 
these basic operations include the Bogoliubov-linear op- 
erations such as the beam-splitter operation, the squeez- 
ing operation and the phase shifter, and nonlinear op- 
erations such as the cross-Kerr interaction ! 36 ! 37 A dis- 
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FIG. 1: Schematic illustration of the dispersive FWM tool- 
box with a dial plate. Two resonator modes Si and a,2 are 
connected to the toolbox. The toolbox is made of a super- 
conducting circuit with controllable parameters ($, $i, <£>2) 
and classical drivings with Rabi frequencies (^1,^2). By ad- 
justing the parameters and drivings to appropriate values, the 
toolbox is switched on to perform selected quantum operation 
on the resonator modes as labelled on the dial plate. 



cussion of these operations can be found in Appendix 
A. In this work, we present a scheme that can gener- 
ate all of these basic operations using one single circuit. 
This toolbox is made of two superconducting qubits cou- 
pled with each other to form a quantum four-level sys- 
tem. Each qubit interacts with one superconducting res- 
onator. By adjusting the parameters of the toolbox, we 
design dispersive four-wave mixing (FWM) processes to 
generate effective quantum operations on the resonator 
modes, as is illustrated in Fig. [TJ During the operations, 
the quantum four-level system is always preserved in its 
quantum ground state by large detunings, and hence this 
scheme is a parametric scheme. Using numerical simula- 
tion, we show that high-fidelity quantum operations can 
be achieved with realistic circuit parameters. In particu- 
lar, we demonstrate the realization of a controlled phase 
gate using the proposed scheme and show that fidelity 
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higher than 99% can be achieved. Compared with previ- 
ous schemesrSr— our proposal provides a switchable cir- 
cuit that can generate all of the basic quantum operations 
by adjusting the parameters of the toolbox. Hence, the 
proposed scheme can advance the scalability of quantum 
information protocols on superconducting resonators by 
connecting a network of resonators with such toolboxes. 

The paper is organized as the following. In Sec. [TTl 
we present the general idea of the toolbox and illustrate 
the idea with a specific circuit made of two coupled su- 
perconducting charge qubits. In Sec. IIIIl we describe the 
dispersive FWM processes for generating effective Hamil- 
tonians on the resonator modes using the toolbox. The 
realizations of four quantum operations are presented in 
detail in Sec. IIVI In Sec. [V] we discuss the main sources 
of quantum errors in this scheme and present our nu- 
merical simulation of the controlled phase gate using the 
effective cross-Kerr interaction. Conclusions are given in 
Sec. IVII In Appendix A, we briefly discuss all of the ba- 
sic quantum operations and their roles in the quantum 
information processing for the resonator modes. 



II. CIRCUIT 

The central element of the toolbox is a quantum four- 
level system that couples with the superconducting res- 
onators and can be constructed in many ways. For the 
discussion in this paper, we will consider a circuit for 
the toolbox that is made of two superconducting qubits 
coupling with each other. The Hamiltonian for the total 
system has the form 



Htot — Hq + H r + H 



(1) 



which includes the Hamiltonian for the four-level system 
= J2j with eigenstates \j) and eigenenergies 

Ej (j — a, b, c, <i), the Hamiltonian for the resonators 
and their couplings to the qubits 




H r = ^ fw ai a}ai + hgi(j xi {a\ + a.,), 



(2) 



FIG. 2: Schematic circuit for the toolbox. The resonators are 
labelled by capacitances CV; and inductances L r i- 



qubit, phase qubit, charge qubit, transmon qubit, and 
etc.J*2- The qubits can be controlled by external electro- 
magnetic fields such as the magnetic flux in the loop of a 
superconducting quantum interference device (SQUID) 
and the bias voltage on a superconducting island, de- 
pending on specific circuit design. Different coupling 
mechanisms between superconducting qubits have also 
been studied, such as capacitive coupling, Josephson cou- 
pling, and inductive coupling^r— Decoherence in super- 
conducting qubits has improved greatly in the past few 
years with T% exceeding 95 /js recently observed42r— 

To illustrate our scheme, we study a toolbox made of 
two superconducting charge qubitsi We want to empha- 
size, however, the FWM approach studied here is a gen- 
eral scheme that can be applied to other superconduct- 
ing qubits with different forms of coupling. As is shown 
in Fig. [21 the charge qubits couple with each other via a 
tunable Josephson junction where the effective Josephson 
energy Ej m can be adjusted by varying the magnetic flux 
$ in the loop. The Josephson energies Eji of the qubit 
junctions can be adjusted by changing the magnetic flux 
$1,2 in the qubit loops. The Hamiltonian for the coupled 
qubits can be derived using a Lagrangian approach with 

H q = (E,n/2)a zl + (E J2 /2)a z2 + H int (3) 
H in t — E mx {a x ia X 2 + b <jyi(jy2 + b a zid z2 ) (4) 



with resonator frequencies u ai and coupling constants 
gi, and the Hamiltonian for the classical drivings on the 
qubits Hp = ^2Kli{t)a X i with driving amplitudes fij. 
Without loss of generality, we have assumed that each 
resonator couples only with one qubit via the a X i opera- 
tor and the classical driving is also on the a X i operator. 
Other forms of coupling and driving could also be con- 
sidered for the toolbox. 

The coupling terms in H r (the driving terms in H p ) 
generate transitions between the eigenstates of the four- 
level system which involve absorption or emission of res- 
onator photons (classical field). The transition matrix 
elements induced by these terms can be derived by pro- 
jecting the a x i operator in the eigenbasis \ j). 

Superconducting qubits in various parameter regimes 
and circuit geometries have been studied, including flux 



where E mx is the charging energy of the capacitance C m 
of the coupling junction and bo = Ej m /AE mx is the ra- 
tio between the Josephson energy and the charging en- 
ergy. Note that the coupling junction is a larger junc- 
tion with its Josephson energy larger than its charging 
energy, in contrast to the junctions of the charge qubits. 
Here, we assume the qubits are biased to have zero charg- 
ing energy so that the qubit energies are the Joseph- 
son energies Eji^ We define the total capacitance con- 
nected to the superconducting island of the ith qubit as 
Csi = Cji+Cgi+C m and the total capacitance connected 
to the ith resonator as C^n = C r i + C g i + Co;, where C g i 
is the capacitance that couples a resonator to its corre- 
sponding qubit, Cji is the Josephson capacitance, C r i 
is the resonator capacitance, and Coi is the capacitance 
that couples to external circuit, as labeled in Fig.[5J With 
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the capacitances satisfying Csri ^> >• C m , the ca- 
pacitive coupling E mx can be derived as 

E mx s» C m e 2 / [CsiCs2 — C m ] ■ (5) 
The eigenenergies of the Hamiltonian H q are 

£ a = E mx b - E s+ /2, (6a) 

E b = -E mx b -E s ^/2, (6b) 

£?o = --E^«6o+-B.-/2, (6c) 

£ d - E mx b + E s+ /2, (6d) 

where £ s± = ^(^jx ± £ J2 ) 2 + 4£^(1 T M 2 can be 
adjusted by varying the qubit energies Eji and the 
Josephson coupling Ej m . The eigenstates are 

| a) =-sin6> + |0 1 2 )+cos6» + |lil 2 ), (7a) 

|6) = cosfl_ Idla) - sin0_ ^Oa) , (7b) 

|c) = sin0_ |0il a ) +cosfl_ IliOa) , (7c) 

|d) = cos0+ |0 x a > +smO+ |lil 2 ) , (7d) 

where |0j) and are the single-qubit eigenstates in the 
cr zi basis and sin 6± = y/[E s ± ± (£./i ± £,/ 2 )] /2£ s ±- 

The resonator-qubit coupling for the ith resonator 
in Eq. ^ can also be derived using the Lagrangian ap- 
proach with 

where C s ^ is the total capacitance of the ith qubit with 
the index i referring to the qubit in the opposite side of 
the toolbox from the ith resonator. We also find that 
the ith resonator couples with the ith qubit as well as 
the ith resonator in the opposite side of the toolbox due 
to the cross-talk of the circuit elements. These indi- 
rect couplings have the forms of hgl cr X i{al + oq) and 

hg^\ai + a{)(d 2 + a 2 ) with coupling constants much 
weaker than the dominant couplings gi in Eq. ([HJ, as 
will be discussed in Sec. [V] 

Both the resonator-qubit coupling and the classical 
driving H p are associated with the qubit operators a X i. 
In the eigenbasis in Eqs. (|7al l7bl I7cl ITdl) . we have 

a xl = cos{9 + - 6-){a ab + a dc ) 

+ sin(0 + - 9-)(<Jdb ~ "WO + h.c. (9a) 
<Jx2 = -sin(6» + + 6-){cr ab - a dc ) 

+ cos(d + + d-)(a ac + (T db ) + h.c. (9b) 

where <xy = |i) (j\ defines the transition operator con- 
necting the eigenstates |i) and \j). From these expres- 
sions, it follows that each resonator (and classical driv- 
ing) generally couples to four transitions. For example, 
the c x i term includes the o abl Odc transitions with am- 
plitude <7i cos(6* + — #_) and the cr ac , Odt, transitions with 
amplitude ±g\ sin(# + — 8 ). 

In the following sections, we will show how to engineer 
the energy levels of the quantum four-level system to have 
the resonators couple only with selected transitions by 
adjusting the circuit parameters. 




bo 

FIG. 3: Eigenenergies of the quantum four-level system with 
circles for E a , up-triangles for E b , squares for E c , and di- 
amonds for Ed- The parameters are Eji/2-Kh = 8.45 GHz, 
Ej2/2nh = 13.95 GHz, and E m x/2itfr = 4 GHz (same param- 
eters as used in the cross-Kerr operation in Sec. II V B|l . 



III. DISPERSIVE FWM SCHEME 

The effective Hamiltonians to implement quantum 
operations on the resonators can be realized via the 
resonator-qubit coupling and classical driving in the 
circuit studied above. Here, we exploit four-photon 
processes^ - — in which single-photon transitions and 
two-photon processes are in the dispersive regime with 
large detunings while the designated four-photon pro- 
cesses are nearly at resonance. As an example, in Fig.@Ji, 
consider a classical driving of frequency ui± generating the 
transition a ac with Rabi frequency f2i and the resonator 
mode at generating the transition adc with effective tran- 
sition matrix element g\ . For the single-photon transition 
induced by the classical driving, the dispersive condition 
requires that |Ai| = \ui — E ca /h\ ^S> fii; and for the 
single- photon transition induced by mode Si, it requires 
that \ui ai — Edc/fi\ 3> 9\\fn\ with Eij = Ei — Ej and 
n\ being the average photon number in a\. For the two 
photon process involving these two transitions, it requires 
that |Ai<5| 3> Qigiy/ni with \S\ = |wi + uj ai — Ed a jh\ be- 
ing the two-photon detuning. Under the dispersive con- 
ditions, the dominant physical processes in this scheme 
are four-photon processes which can generate effective 
coupling between the resonator modes. Because the dis- 
persive conditions prevent real transitions, the quantum 
four-level system is preserved in its ground state during 
the operation. The processes studied here are hence para- 
metric schemes where the four-level system is subject to 
a "quantum" energy shift 4^ 

To implement a quantum operation, we need to adjust 
the parameters of the superconducting circuit to find ap- 
propriate effective coupling constants and energy separa- 
tions between the eigenstates. For the circuit in Fig. [2] 
we adopt two approaches to determine the circuit param- 
eters. The first approach uses the relation in Eq. (I9"al 
I9bl) to tune the effective transition matrix elements by 
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adjusting the angles 9±. By adjusting the circuit pa- 
rameters Eji , Ej2 , b , the angles can be varied in a large 
range. The second approach exploits the controllability 
of the energy levels to engineer large detunings to sup- 
press unwanted transitions. As an example, in Fig. [3l 
we plot the eigenenergies as functions of the ratio bo (de- 
fined in Sec. |TTJ) . As bo varies, the state \a) (|c)) can be 
either above or below the state \b) (\d)), which provides 
the possibility to arrange both the order and the energy 
separation of the eigenstates. To make, e.g. the mode 
hi only couples strongly to the a ab transition, we choose 
the parameters to have: (1) cos(6> + — 6 ) 3> sin(6'+ — 9 ) 
and (2) \uj ai — E ba /h\ <C |w 0l — E dc /h\. The first condi- 
tion significantly reduces the effective coupling between 
di and the <J ac ,Odb transitions and the second condition 
suppresses the h\Udc term by large detuning. Now, mode 
di couples mainly with the a a b transition. In general, by 
combining the above two approaches, all quantum oper- 
ations described in Appendix A can be realized with the 
FWM scheme within one single circuit. 



IV. REALIZATION OF OPERATIONS 

In this section, we present the realization of four quan- 
tum operations described in Appendix A using the quan- 
tum toolbox. The system parameters are chosen as the 
following: ui ai /2ir = 10 GHz and uj a2 /2ir = 16 GHz for 
the resonator frequencies, gif2.ii = gijlix = 0.3 GHz 
for the coupling constants, and E mx /2irh = 4 GHz for 
the capacitive coupling defined in Eq. (JSJ) - The effective 
Josephson energies Ej%, Ej2, and Ej m can be adjusted 
by tuning the magnetic flux $i, $2, and $ respectively. 
The energy diagrams of the dispersive FWM schemes for 
the operations are shown in Fig. 2J 



A. Beam-splitter operation 

Consider the energy diagram in Fig. 0^, where the de- 
tunings are defined as Ai = uii — E ca /h, A2 = uj2—E ba /h, 
and S = u>i+ u> ai — Eda/h with being the frequency of 
the classical driving. In the energy diagram, the res- 
onator modes only couple with the upper transitions 
Cdb, Ode and the classical drivings only couple with the 
lower transitions a abl cr ac . This can be achieved by choos- 
ing the parameters bo = —0.68, Eji/2irh = 8 GHz, and 
Ej2/2irh =15 GHz so that the large separation between 
the ground state and the excited states provides energy 
selection for the transitions. Under these parameters and 
appropriate classical driving frequencies, the detunings 
for the desired single-photon and two-photon processes 
as labelled in the energy diagram are A,/27r = —4 GHz 
and 8/2tt — 3.49 GHz. The detunings for the unwanted 
transitions, e.g. U) ax — E ba /h,uj a2 — E ca /h,~ 9 GHz, are 
much larger so that the unwanted transitions are sup- 
pressed. These parameters also give | cos(#+ — ^> 
I sin(# + — 6-)\ so that the coupling between the mode 




FIG. 4: (Color online) Energy diagram for the FWM schemes 
of (a) beam-splitter operation, (b) cross-Kerr interaction, (c) 
two-mode squeezing, and (d) single-mode squeezing. The la- 
bels are: detuning Ai for single-photon transition, detuning 
5 for two-photon process, detuning Af for four-photon pro- 
cess, Rabi frequency Q,i for classical driving (red arrows), and 
operator ca for resonator mode (light and dark blue arrows). 

di (0,2) and the transition tjdb (ede) is much weaker than 
the coupling between the mode di (02) and the transi- 
tion a dc {(Jdb)- As a result, the mode di mainly couples 
with the adc transition. Similar arguments apply to other 
transitions in the energy diagram. 

We divide the total Hamiltonian into two parts: H to t — 
Ho + V where 

Ho/H = OJi <T CC +LJ 2 &bb + {Wl +U ai ) Odd +^ ^ai «1 «i ) (10) 

i 

f« — an( i ^ includes all remaining terms in the 

total Hamiltonian. Here, the ground state energy is set 
to zero. The Hamiltonian in the interaction picture of Ho 
can be written as e iH t/hy e -iH t/h ^ Hr) + Vj UIK j e r the 

rotating wave approximation with Hio/h = — Ai<j cc — 
A 2 cr 6b - 5a d d and 

Vj/H = Vli(j ca + Vt 2 (Tba + giaiVdc + g2e tAFt a 2 (Tdb + h.c. 

(11) 

The effective couplings gi — gicos(0 + — 6*_) and 02 = 
g2COs(9 + + 0-) are derived from Eqs. (|9"al I9b[) . And, 
Ap = lji + ui ai — L02 — uj a2 is a small detuning for the 
four-photon process designed to balance the extra terms 
in the effective Hamiltonian for the resonators. 

Given the dispersive conditions discussed in Sec. lIIIi we 
treat Vi as a perturbation to the Hamiltonian Hjo- As- 
sume that the toolbox is initially prepared in its ground 
state I a). It can be shown that the dominant correction 
to Hjo by the perturbation Vi is a fourth-order term that 
generates a "quantum" energy shift in the ground state 
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in the form of (7 aa H^ m where 
H hm /h = 5e\ m a\ai + (x hm a\a 2 e 1 ^ 1 + h.c)j (12) 

i 

with the energy shifts Se° m = Qfgf/AfS, the effective 
coupling constant 

x hm = fiifi 2 5ig 2 /A 1 A 2( 5, (13) 

and the four-photon detuning Ap = 8e\ m — Se\ m . The 
small four-photon detuning is chosen to balance the ef- 
fect of the mode shifts 8t° m . The above effective Hamil- 
tonian performs the beam-splitter operation on the res- 
onators while the toolbox is preserved in the ground state 
during the operation. The beam-splitter operation can 
perform a swap gate on the resonators after applied for 
a gate time t g — 7r/2|x bm |'^ With the parameters given 
above and with Sli/2n = 1.5 GHz, the effective coupling 
is |x bm |/2"" = 6.2 MHz and the swap gate has a gate time 
of 40.5 ns. 



B. Cross-Kerr nonlinearity 

The energy diagram for the cross-Kerr operation is 
shown in Fig. |4Jd with the detunings defined as Ai = 
uJ ai -E ba /h, A 2 = uj a2 -E ca /h, and S = u> ai +uj a2 -E da /h. 
Compared with the beam-splitter operation, the asym- 
metry in the energy levels is reduced and no classi- 
cal driving needs to be applied. The parameters are 
chosen to be bo = —0.61, Eji/2nh = 8.45 GHz, and 
E J2 /2Trh = 13.95 GHz, which give Ai/2tt = -4.59 GHz, 
A 2 /2tt = -4.93 GHz, and 6/2n = 0.17 GHz. Under these 
parameters, the effective coupling constants for the de- 
sired couplings aiCTdc, aicr a (, are gi/2n « 0.3 GHz. While 
the unwanted couplings aiadb, aic QC are strongly sup- 
pressed with <?i/27r w 0.01 GHz. The unwanted coupling 
a 2 Od c has a detuning of \ui a2 — Ed c /fi\/2iT ss 11 GHz and 
is suppressed by the large detuning. 

We divide the total Hamiltonian as H to t = Ho + V 
with 

H /h = UJ ai a bb + LO a2 ^cc + U t Odd + '^Lo . l a\a l (14) 

i 

and uJt = 0J ai + uj a2 . The Hamiltonian in the interaction 
picture of H is e *H t/hy e -iH t/h ra Hiq + Vj UIK } er the 

rotating wave approximation where Hiq/H = — Aiay, — 

A 2 er cc - 5<7dd and 

Vj/H= g\ai(<7 dc . + <7ba) + h^i^db + <7c.a) + h.c. (15) 

with the effective couplings g\ = g\ cos(#+ — 0_) and 
92 = 92 cos(8 + +8-). With the same perturbation theory 
approach as used in the previous subsection, we derive 
the fourth-order perturbation correction to the ground 
state energy a aa H c ^. Here, 

#ck A = 6e t4*i + X ck a\4a2ai , (16) 




£(ns) 

FIG. 5: Time evolution of \{ipi\ip{t))\ under the cross-Kerr 
interaction. Solid curve: is the target state at time 

t g \ dashed curve: \ipi) is the initial state given in Eq. (|19[1 ; 
thin-solid curve: \ipi) = |a) IO1O2) ; and dotted curve: = 
|o>|lil a >. 

with the energy shifts Sef* = fff/Aj and the effective 
coupling 

x ck = (1/Aa + 1/A 2 ) 2 {gtgl/8) . (17) 

With the above parameters, |x ck |/27r = 6.3 MHz. 

Given the effective cross-Kerr interaction in Eq. (|16l) . 
the controlled phase gate can be realized on the resonator 
modes which generates the following transformations 

|0i0 3 ) -> |0i0 2 ), (18a) 

|0il 2 )^e^ kt |0 1 l 2 ), (18b) 

|li0 3 > -> e-^lliOa), (18c) 

|lil a ) e - 4,5e i kt -^ kt - lxCkt |lil 2 ), (18d) 

on the corresponding resonator states. At t g = 7r/|x ck |i 
the nonlinear term generates a phase n on the state 
1 1 1 1 2) - With the above parameters, we find that the 
gate time is t g = 79.4 ns. To test the above results, we 
simulate the time evolution of the combined resonator- 
toolbox system under the Hamiltonian Hiq + Vj given in 
Eq. (fT5")) . Consider the initial state 

|a)(|0 1 > + |l 1 »(|0 2 ) + |l 2 ))/2, (19) 

where the toolbox is in the ground state \a). The target 
state at time t g can be derived by substituting time t with 
t g in Eqs. (|18a|18bll8cll8dj) while keeping the toolb ox in 
the state \a). In Fig. [51 we plot \(ipi\ip(t))\ for different 
states \ipi). After the evolution starts, the amplitude of 
the initial state starts to decrease. At the target time t g , 
the amplitude of the target state nearly reaches unity, 
which demonstrates that the final state of the system 
has successfully evolved to become the target state. The 
amplitude of the state |a)|0i0 2 ) stays at 0.5 with no sig- 
nificant transition to other states during the evolution. 
Meanwhile, the amplitude of the state |a)|lil 2 ) shows 
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small deviation from 0.5 due to leakage to the excited 
states of the toolbox during the evolution, which can re- 
duce the fidelity of the controlled phase gate. 



C. Two- mode squeezing 

In Fig. 0J:, we present an energy diagram that can real- 
ize the two-mode squeezing operation with the detunings 
defined as Ai = uii — E a b/h, A 2 = 102 — Edb/h~, and 
S = u>2 — co ai — E c b/h. Here, the order of the eigen- 
states is switched with E a > Ef,, i.e. the state \b) is 
now the ground state of the toolbox. The circuit param- 
eters are chosen to be bo = 1.2, Eji/2ttH = 8 GHz and 
Ej2/2ttH — 14 GHz. Under these parameters, we tune 
the classical driving frequencies to give Ai/27r = 3 GHz, 
A 2 /2tt = —5 GHz and 6/2ir = -3.67 GHz. The clas- 
sical drivings induce the <J a b,<Jdb transitions in the en- 
ergy diagram. It can also be shown that the mode a± 
mainly couples with the <idc transition and the mode 62 
mainly couples with the a ac transition. With these pa- 
rameters, the driving frequencies are uji/2ti = 10.9 GHz 
and W2/27T ~ 24.9 GHz. Note that at high driving fre- 
quency, the microwave driving can induce quasiparticle 
excitations in the superconductor. This effect can be 
avoided by considering a two-photon process for the <Jdb 
transition. 

Using the approach in the previous subsections with 

H Q /H = LOiCTaa + U S (J CC + L0 2 ^dd +^^4"'' ( 20 ) 

i 

and uj s = ll>2 — w ail we derive that Hio/h — — Ai<r aa — 
A 2 cr dd - 5a cc and 

Vj/h = Qi<7 ab + Q2&db + 9\a\<J dc + 32e lAFt a 2 cr ac + h.c. 

(21) 

with the effective couplings <?i = g\ cos (9 + — 6U) and 
.92 = 92 cos(0+ + 0-), and a small four-photon detuning 
Ap = uj2 — uj\— cl> Qi — u)a 2 ■ The fourth-order perturbation 
correction to the ground state energy is cr aa H sq with 

H sq /h = Sd?alai + (x sq aI4e lAFt + h.c.') , (22) 

i 

where the energy shifts are <5e^ q = il?gf/A?S with the 
index i referring to the circuit elements in the opposite 
side of the ith resonator (e.g. for i = 1, i = 2), the 
effective coupling is 

x sq = QiQaSi^/AxA^, (23) 

and the four-photon detuning is Ap = — Se^ 1 — 5e s 2 q . 
The nonzero four-photon detuning balances the effect 
of the energy shifts and makes the squeezing operation 
possible. With the parameters given above and with 
QipTr = 2 GHz, we have |x sq l/27r = 4.3 MHz. 



D. Single-mode squeezing 

The single-mode squeezing operation can be realized 
with the energy diagram in Fig. |4ji with the detunings 
defined as Ai = ui ai — E a b/h, A 2 = u>2 — Edb/h, and 
S = UJ2 — w ai — E c b I h. The circuit parameters are cho- 
sen to be & = 1-2, Ejx/2-kH = 6 GHz and Ej 2 /2Tvh = 
15 GHz. Tuning the frequencies of the classical driv- 
ings, we have Ai/2-k = 3 GHz, A2/27T = —5 GHz, and 
5/2ir = —4.75 GHz. The circuit parameters are chosen 
to have the energy separation E a b close to the frequency 
of mode a\. Hence, a\ couples strongly with the a a b 
transition as well as the a dc transition. Meanwhile, the 
energy levels are adjusted so that the frequency of mode 
a,2 is largely detuned from all possible transitions and 
is effectively decoupled from the toolbox. The classical 
drivings generate the cr ac , Odb transitions. With the above 
parameters, ui/2ir = 12.0 GHz and W2/27T = 25.0 GHz. 

With the Hamiltonian 

H /h = UJ ai (Jaa +uj s cr cc + uj2^dd + ^ uj ai a\ai, (24) 

i 

and uj s = uj 2 - w ax , we have H IQ /h = -Aia aa - A 2 a d d - 
Sa cc and 

Vi/h = n 1 a ca e lAFt + n2<r d b + giai(crdc + crab) + h.c. (25) 

with the effective coupling constant g± — g\ cos(0 + — 6 J) 
and the four photon detuning Ap = LU2 — u)\ — 2oj ai . 
Then, we derive the fourth-order perturbation correction 
to the ground state energy a aa H s(l i with 

H sql /h = Se sql a\ai + (x sql a\a\e- iAFt + h.c?) (26) 
where the energy shifts are 

^ ql - (8/A 1 + n?/A? + nj/Aj) (gl/S) , (27) 
the effective coupling constant is 

X^ 1 = il^gl/A^S, (28) 

and the four-photon detuning is Ap = 2Se bql . With the 
parameters given above and with D,i/2n — 1 GHz, we 
have |x sql |/27r = 11.1 MHz. 

V. ERROR SOURCES 

Quantum errors can affect the effective quantum oper- 
ations on the superconducting resonators. In our system, 
the main sources of quantum errors include (1) unwanted 
transitions induced by resonator-qubit coupling and clas- 
sical driving, (2) indirect coupling due to the cross-talk 
between different circuit elements, and (3) decoherence 
of the qubits and resonators. 

We first study the effect of unwanted transitions on 
the effective quantum operations studied in Sec. IIVI In 
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7 (3) 



(ai + a\)(a,2 + a\) with the coupling constant 



E mx /2ir (GHz) 



FIG. 6: Fidelity and gate time versus E mx for the controlled 
phase gate using parameters in Sec. IIVBI 



the proposed scheme, the circuit parameters can be ad- 
justed to suppress the unwanted transitions either by 
reducing the coupling matrix element of the unwanted 
transitions or by varying the energy levels of the un- 
wanted transitions to produce large detuning. To test 
the effectiveness of our approach, we numerically simu- 
late the quantum operations using the full Hamiltonian 
in Eq. ([1]) which includes all of the nonzero transition 
matrix elements. In Fig. [Bl we present our result for 
the fidelity of the controlled phase gate which is demon- 
strated in Sec. IIVBI The initial state for the simulation 
is given in Eq. (|19p . For each value of the coupling en- 
ergy E mx with E mx /2n € (3.5, 4.5) GHz, we search for 
the maximum fidelity by varying Eji, bo, and the gate 
time. The maximum fidelity and the corresponding gate 
time are then plotted in Fig. [5] We show that the max- 
imum fidelity can exceed 0.99 at E mx /2n — 4 GHz. Our 
simulation hence shows that high fidelity can be achieved 
for effective quantum operations by choosing appropriate 
parameters. In Fig. [5J the fidelity fluctuates "randomly" 
within a narrow range as E mx varies. This is because the 
time evolution of the total system includes small but fast 
oscillations resulted from multiple unwanted off-resonant 
transitions, as can be seen in Fig. [3] Such oscillations 
affect the optimal gate time and the fidelity by a small 
magnitude in a nearly "random" way. 

Another source of quantum errors is the circuit cross- 
talk between different circuit elements. In Eqs. (JUIH]), the 
direct coupling between the resonator and its neighbor- 
ing qubit is presented. However, due to the cross-talk, 
each resonator also couples with the qubit in the opposite 

side of the circuit in the form of %j 2 a x i(al + aj). The 
coupling constant for this indirect coupling is 



-,(2) 



CeiCs2 — 



2C- 



-hljj n 



En 



1/2 



(29) 
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(3) _ \/ Cg\Cg2C m I Cg\Cg2 



CeiCe2 ~ V 4CE r iCsr2 

It can be shown that 



yj tvjJ ai h0J a2 ■ (30) 



Cm 



C a 



tlWn 



4C* S „ V e 2 /2C g 



< 1. 



(31) 



These indirect coupling terms can hence be neglected. 

Decoherence of the superconducting qubits and res- 
onators is one of the key barriers for scalable quantum 
information processing and has been intensively studied 
in the past two decadesi 44 i 46 The quantum toolbox stud- 
ied here is made of two coupled superconducting qubits 
and can be subject to both single-qubit and two-qubit 
decoherence. However, for the dispersive FWM scheme 
proposed in this work, the quantum toolbox is largely 
preserved in its ground state during the quantum opera- 
tions, and hence, is not affected by either single-qubit or 
two-qubit decoherence. Only small leakage to the excited 
states can be induced by the unwanted transitions as is 
studied above. Given the small amplitude of the leakage 
as is shown in Fig. [5j qubit decoherence will not affect 
the effective quantum operations significantly. Super- 
conducting resonators can have relatively high Q-factors. 
With Q exceeding 10 6 , it corresponds to a damping time 
of ~ 100 fis. Meanwhile, our study in the previous section 
shows that the time scale for the quantum operations - 
~ ^/Ix"! f° r operation a - is below 100 ns which is shorter 
than the decoherence time by 2 — 3 orders of magnitude. 



VI. CONCLUSIONS 

To conclude, we presented a superconducting quantum 
toolbox that can perform various quantum operations on 
superconducting resonators in one single circuit. The 
scheme exploits the dispersive FWM approach to gen- 
erate effective couplings between the resonator modes. 
By adjusting the circuit parameters, the energy levels 
and coupling constants can be varied to generate spe- 
cific quantum operation. We discussed the main error 
sources in the schemes and numerically simulated the 
controlled phase gate on Fock states. Our results showed 
that high-fidelity quantum operations can be achieved 
in this circuit. One advantage of this scheme is that 
nearly all quantum operations for both discrete-state and 
continuous- variable quantum protocols can be realized in 
one single circuiti 36 i 37 ' 54 Our scheme can advance the im- 
plementation of quantum information processing on the 
microwave modes in the superconducting circuits. 
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APPENDIX A 

Quantum information processing with photons has 
been widely studied either with discrete states or con- 
tinuous variable states . 36 ' 3 Here, we briefly summarize 
the basic quantum operations for the photon modest 
Two categories of quantum operations are considered: 
the Bogoliubov-linear operations and the nonlinear in- 
teractions. The photon modes are represented by the 
annihilation (creation) operators a\ and 02 (a\ and d\). 

A. Bogoliubov-linear operations 

The Bogoliubov-linear operations perform the follow- 
ing transformation 

o< -> + B a h ) + c ' ; ( 32 ) 

3 

with coefficients Aij , Bij d . An arbitrary Bogoliubov- 
linear operation can be constructed using the basic ele- 
ments: the beam-splitter operation, the squeezing oper- 
ation, and the phase shifter. 

The beam-splitter operation can be realized by i?bm = 
JiX hm e l ^a\a2 + h.c. with coupling amplitude x° m . Under 
this Hamiltonian, the operators evolve as 

\a 2 (t) J ~\ -e-^sixnp cos ip ) \ a 2 (0) ) 

(33) 



with (p = Xbm^- At ip = 7r/2, the beam-splitter operation 
swaps the states of the two modes up to a phase fac- 
tor. In discrete-state quantum computing schemes, this 
operation can generate single-qubit Hadamard gate. 

The squeezing operation can be realized by H sq = 
ihx sq a\a 2 + h.c. with coupling amplitude x sq . Under this 
Hamiltonian, the operators evolve as 

( ai{£) \ _ ( cosh ip sinh^ \ / a x (0) \ 

\a\it) ) \ sinht^ cosht^ J \at{0) J { ' 

with ip = x sq ^: which describes the parametric amplifica- 
tion process that generates two-mode squeezing.-^ When 
applied to the vacuum state, it generates the so-called 
continuous-variable EPR states. When combined with 
the beam-splitter operation, it can generate squeezing 
on individual mode. A related operation is the single- 
mode squeezing operation which can be generated by 
H sq i — ihx sql {a^)1 + h.c. on mode a,-. 

The phase shifter operation can be realized by H p h = 
hA ph a\ai which creates a shift in the resonator frequency. 
The above quantum operations can be combined to gen- 
erate arbitrary linear transformations in Eq. (|32[) . 
B. Cross-Kerr nonlinearity 

One nonlinear operation is the cross-Kerr interaction 
given by H c ^ — hx ck a\did\a,2 between two modes with 
interaction amplitude x ck - This interaction can lead 
to controlled gates on photon qubits.— For continuous- 
variable schemes, this operation together with the lin- 
ear operations can generate operations that are arbitrary 
polynomials of the quadrature variables i 36 ' 37 ' 48 This op- 
eration can also be exploited for quantum nondemolition 
measurement on photon states. 
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